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Motivation 


•  How  will  energy  extraction/introduction  via  joule 
heating(J  •  E)  affect  the  stability  of  the  detonation? 

•  How  will  the  magnetic  Pressure(J  x  B)  evolution  affect  the 
Reaction  Zone  and  detonation  stability? 
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Goals  of  the  Present  Study 


•  In  order  to  understand  detonation-magnetic  field  interactions, 
one  must  first  understand  the  stability  criteria  of  an 
unmagnetized,  unsupported  detonation. 

•  Explore  the  nonlinear  dynamics  involved  in  detonation  stability 

►  Induction  lengths  relation  to  kinetics  &  dynamics 

•  Examine  the  coupling  of  large  and  small  length  scale  physics 

►  Correlating  different  modes  of  peak  pressure  behavior  to  small 
scale  phenomena 
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Numerical  Methodology 


Inviscid,  one-dimensional  Euler  equations  using 
multi-step, reversible  reaction  mechanism: 

Qt  +  F(Q)X  =  S(Q)  (1) 

where  the  vectors  represented  by  Q,  F,  and  S  are,  respectively, 


Q=  (  pi  V  F=  (  pu2  +  p  \  ,  S 

\  E  J  V  (E  +  P)u  ) 


Us  \ 

0  (2) 
Es^eos  J 


where  the  total  mixture  density  p  =  J2sPs'  ar|d  the  total  energy  E 
may  be  written 


cv{  T)dT  +  -pu2 


(3) 


Distribution  A:  For  Public  Release;  Distribution  Unlimited 


Goals  Validation  Simple  Det.  Grid  Conv.  Ind.  Zone  Model  Euler  Ideal  MHD  Kinetics 

Numerical  Schemes 


•  Monoticity  Preserving(MP)  Schemes  (Suresh  &  Huynh,  1997) 

►  5th  order  spatial  discretization  was  used  in  conjuction  with  3rd 
order  TVD-Runge-Kutta  time  integration 

►  Contact  Discontinuities  well  resolved  without  the  use  of 
artificial  compression  methods 

•  Advection-Diffusion-Reaction  Weighted  Essential 
Non-Oscillatory(ADERWENO)  Scheme  (Titarev  &  Toro, 

2001) 

►  5th  order  spatial  and  3rd  order  temporal  without  Runge-Kutta 
time  integration 

►  Utilizes  Lax-Wendroff  procedure  and  Taylor  series  expansion  of 
WENO  fluxes  for  high  order  in  time 
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•  Euler  equation  solutions  using  MP  schemes  validated  for 
standard  problems: 

►  Sod’s  ID  Shock  tube  problem 

►  Lax's  ID  problem,  shock  tube  with  velocity  field 

►  Shu-Osher  ID  problem,  entropy  wave-shock  interaction 

►  ID  Blastwave  problem,  e.g., 


-  Exact 


•  MP5 


300  grid  points 


X  [m\ 
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Ideal  MHD  Conservation  Equations 


Q  t  +  Fx  —  0 


/  Ps  \ 

/  PsUn  \ 

pu 

,  F  = 

puun  +  P*n  —  B6n 

B 

Ur B  -  uBn 

V  E  J 

\  (E  +  P*)  un  —  Bn(u  ■  B)  / 

where 

Bn  =  Bx  nx  +  By  riy  +  Bz  nz,  P  =  P  +  ^||Bj| 

Validated  for  standard  ID  Brio-Wu  Problem,  analogous  to  Sod’s 
shock  tube: 


PL,  Pl,  ul 


B, 


-By  Pr,Pr,ur 
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ID  Brio-Wu  validation  with  ADERWENO 


x  M 


Initial  Conditions: 

{pl,Pl,ul,Bx,By} 

{PR,PR,UR,BX,By} 


=  {1^  ,  105  Pa  ,  Of  ,0.75,1} 

=  {0.125^  ,  104  Pa  ,  Of  ,0.75, -1} 
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Chemical  Kinetics  Conservation  Equations 


Operator-Splitting 

dQ=s 

dt 

where 

(  Ps  ' 

\  /  OJs 

Q  = 

pu 

•  s  =  0 

! 

\  E  , 

/  \  E^seos 

Ws  =  ^  Vrskfr  ^  Urskbr  JJjXjf'-/ 

r  j  r  j 

"rk  =  v"k  -  v'rk 

v"k\  coefficient  of  kth  species  in  the  rth  forward  reaction 
v'rk\  coefficient  of  kth  species  in  the  rth  reverse  reaction 
[Xs]:  Concentration  of  sth  species 
eos:  Internal  Energy  of  formation  at  OK 
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Kinetics  of  Hydrogen-Air  Mixture: 

•  The  chemistry  includes  eight  reacting  species,  /-/2,  O2,  H,  O, 
OH,  HO2,  H2O2,  H2O,  and  the  non-reacting  diluent  A/2. 

•  Thirty  eight  elementary  reactions  are  used  in  this  mechanism 
and  the  backward  rates  are  computed  from  equilibrium 
constants. 

•  Convection  and  Kinetics  were  operator  split 

►  Point  Implicit  Euler  was  used  to  solve  for  the  kinetics 
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Induction  Delay  Time 


Reduced  H2  —  O2  —  A/2  Reaction  Kinetics  (9  species,  38  reactions) 


Can  be  fitted 
as  t  ~  a(P)e  r 


Criterion  for  Induction:  Maximum  atomic  Hydrogen  concentration 
Mixture:  Stochiometric  H2  —  O2  —  A/2 
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Detonation  Test  Setup 


P  =  40  atm 
T  =  1500  K 


L 

K - 30  cm - * 

•  Premixed  Stoichiometric  Mixture  of  AV2— Air 

•  Closed  Ends 

•  Spark  ignited  (L  =  0.25  cm  in  the  Present  Study) 

•  Dcj  ~  2054  m/s 


P  =  1  atm 
T  =  300K 
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Spark  Ignited  Detonation  -  Pressure  Contour  (MP5) 
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Spark  Alteration  ( Pspark  =  50  atm,  L  =  ,25cm),  MP5 


1  ^  r3  4  Xfcml  6  7  8  r 

Grid  Resolution(Ax)  =  5/xm  Dcj  =  2054™ 

L  =  0.25  cm  D  =  203l| 
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Typical  Peak  Pressure  vs  Time  Plot 


Grid  Resolution:  Ax  =  2.5 [im 

HF-  High  Frequency,  HA-  High  Amplitude 
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Time  to  Re-Explosion 


Time  of  re-explosion (Texp): 
•  criteria  P  >  100  atm 
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Grid  Convergence 


•  Peak  Pressure  vs  Time  data  relates  macroscopic  phenomena 
to  microscopic  phenomena 

•  Spectral  content  of  the  High  Frequency  and  High  Amplitude 
Modes  for  various  grid  resolution  can  be  used  to  determine 
convergence 

•  High  Frequency  modes  were  inconsistent  as  the  grid  resolution 
increased  to  Ax  >  7.5 /im 
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Typical  Peak  Pressure  vs  Time  Plot 


Grid  Resolution:  Ax  =  12.5 fim 
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Spectral  Convergence 
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Simplified  Model 

Induction  Zone  Dynamics 


flame 


Induction  Zone 

*  _ 


shock 


single  period 


£ 

dx 

dt 


dx 

dt 


=  tv2(x,  t) 

entropy 

=  c(x,  t )  -  u2(x,  t) 


acoustic 

where  U2(x,  t)  =  |u(x,  t)  —  D(t)|  (detonation  ref.  frame) 


Entropy 

Acoustic 


Wave 


(4) 

(5) 
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Verify  Simple  Model 


f  =  T  1 ,  frequency 
Assumptions  for  zeroth  order  approximation: 

.  D  =  y  $  D(t)dt 


•  7  ~  1.28 

•  T,  P.  p,  u  — >  F(x,  t), 


9F{x,t)  ~ 
dt 


o  & 


~o 

ox 


(6) 


Distribution  A:  For  Public  Release;  Distribution  Unlimited 


Temperature  [K] 


Goals  Validation  Simple  Det.  Grid  Conv.  Ind.  Zone  Model 


Verify  Simple  Model 

High  Frequency  Mode 


a:  29.1/iS  b:  29.3 fis  c:  29.5 /xs 
f  ~  2.08  Mhz,  In  good  agreement  with  spectral  analysis 
(2.29  ±  0.4  Mhz) 
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Verify  Simple  Model 

High  Amplitude  Mode 


a:  56.7 /is  b:  58.3 /is  c:  60.1/rs 

f  ~  310  khz,  In  great  agreement  with  spectral  analysis 
(310  ±  40  khz) 
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Same  fundamental  dynamics  for  High  Frequency  and  High 
Amplitude  Modes 

The  location  of  the  hot  spot,  whether  within  the  flame  or 
reaction  zone,  plays  a  key  role  in  pressure  oscillations. 

►  ‘hot  spot'  inhibits  progress  of  flame  toward  shock  by 
pre-igniting  fluid,  thus  suppressing  the  peak  amplitude  of 
pressure 

HF  Mode:  “hotspot"  well  resolved  within  induction  zone, 
leading  to  small  fluctuations  in  Pressure 

HA  Mode:  fluctuations  are  within  flame,  allowing  for  more 
energy  release  (via  Swacer  effect) 
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Future  Work 


•  Conduct  similar  test  with  a  seeded  species  of  low  ionization 
energy  and  direct  initiation  (Hydrocarbon  fuel) 

►  Observe  how  ionization  processes  effect  the  induction  region 
and  large  scale  phenomena 

•  Apply  B-field  of  varying  strength 
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Backup 

MHD  Effects 


Additional  Physics  added  in  MHD 

F  =  J  x  B  Lorentz  Force 

_  (7 

Q  =  J  ■  E  Joule  Heating 

Questions: 

•  How  will  energy  extraction/introduction  via  joule  heating 
effect  the  stability  of  the  detonation? 

•  How  will  the  magnetic  Pressure  evolution  effect  the  Reaction 
Zone  &  detonation  stability? 
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